Draft version February 1, 2008 

Preprint typeset using IAT^X style emulateapj v. 04/03/99 



o 
o 
o 

(N 

Or 

<D 

t> 

(N 
> 

00 
00 

o 
o> 
o 
o 

O 



Oh: 

6 

C3 



X 



TWO-DIMENSIONAL MHD NUMERICAL SIMULATIONS OF 

MAGNETIC RECONNECTION TRIGGERED BY 

A SUPERNOVA SHOCK IN INTERSTELLAR MEDIUM: 

GENERATION OF X-RAY GAS IN GALAXY 

SYUNITI TANUMA 1 ' 2 ' 3 , Takaaki Yokoyama, 4 , Takahiro Kudoh 2 ' 5 , 

AND 

Kazunari Shibata 6 
Draft version February 1, 2008 

ABSTRACT 

We examine the magnetic reconnection triggered by a supernova (or a point explosion) in 
interstellar medium, by performing two-dimensional resistive magnetohydrodynamic (MHD) 
numerical simulations with high spatial resolution. We found that the magnetic reconnection 
starts long after a supernova shock (fast-mode MHD shock) passes a current sheet. The current 
sheet evolves as follows: (i) Tearing-mode instability is excited by the supernova shock, and the 
current sheet becomes thin in its nonlinear stage, (ii) The current-sheet thinning is saturated 
when the current-sheet thickness becomes comparable to that of Sweet-Parker current sheet. 
After that, Sweet-Parker type reconnection starts, and the current-sheet length increases, (hi) 
"Secondary tearing- mode instability" occurs in the thin Sweet-Parker current sheet, (iv) As 
a result, further current-sheet thinning occurs and anomalous resistivity sets in, because gas 
density decreases in the current sheet. Petschek type reconnection starts and heats interstellar 
gas. Magnetic energy is released quickly while magnetic islands are moving in the current sheet 
during Petschek type reconnection. The released magnetic energy is determined by the inter- 
stellar magnetic held strength, not energy of initial explosion nor distance to explosion. We 
suggest that magnetic reconnection is a possible mechanism to generate X-ray gas in Galaxy. 



Subject headings: magnetic reconnection — Galaxy: structure 
MHD — shock waves — X-rays: ISM 



ISM: magnetic helds 



1. INTRODUCTION 

Observations of the diffuse component of X-rays 
(e.g. Snowden 1995) indicate that hot components of 
temperatures higher than 10 6 K occupy a consider- 
able volume of the interstellar space. These hot com- 
ponents may be generated by supernova remnants 
(SNRs), stellar winds and magnetic heating. X-rays 
. from Galactic ridge (GRXE; Galactic Ridge X^ra> 



space in this paper. The mean strength of interstel- 
lar magnetic field i s several /iG (|Beck et al. 1996 ; 
Heiles et al.~1976t |Sofue et al. 1986; ; |Vallee 19971) . 



The interstellar magnetic field has random compo- 
nents and generates possibly many current sheets at 
all scales (e.g. |Ohno k, Shibata 1993Q . Parker (1992) 
pointed out the importance of magnetic reconnection 
for the heating of Galactic plasma (see also Hanasz &| 



1^ ° ""' ■ i ■' " ° ■ v frr ' i 7 ",~ -i na^lr^ 5 Lesch 1998J) . Magnetic reconnection is a fundamental 
Emission; see discussion; Kanoda et al. 1997 j- Koyama ; : — : — -i' ° . 



et al. 1986a|; |Makishima 1994|; |Yamasaki et al. 1997| ; 
Yamauchi et al. 1996| ; |Valinia fc Marshall 1998| ), for 
example, can not be fully explained by the point 



sourc es such as SNRs and stellar winds ( Kanedaet al. 
1997| ) . We study the magnetic heating in interstellar 



intrinsic property of agitated magnetized, turbulent 
plasma. Whenever the magnetic fields collide with 
another field with different direction by, for example, 
the shear motion, magnetic energy is rapidly dissi- 
pated. Outside the forming current sheet the mag- 
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netic field lines are frozen- in to the interstellar gas, 
which is provided by the high electrical conductiv- 
ity. When the oppositely directed field lines collide 
with each other, the field gradient steepens and the 
current density increases until strong dissipation sets 
in (e.g. via anomalous resistivity) to trigger fast re- 
connection (e.g. Hanasz & Lesch 1998; Ugai 1986J ). 
The reconnection heats the interstellar gas by releas- 
ing interstellar magnetic energy, and accelerates it by 



forming two dimensional (2D) numerical magnetohy- 
drodynamic (MHD) simulations including the elec- 
tric resistivity. Many MHD numerical simulations 
have been performed for magnetic reconnection in the 
solar atmosphere (e.g. Yokoyama & Shibata 1995; 



Odstrcil fc Karlicky 1997 ). Some similar attempts 
have been made in the model of intracluster medium 
([Valinia et al. 199€; see also|Tajima fc Shibata 1997) 



mag netic tension force to Alfven velocity, forming in- 



terst ellar jets. The magnetic reconnection plays an tostar ( Hayashi, Shibata, & Matsumoto 199£; Shi- 



Galactic halo (Birk, Lesch, & Neukirch 1995; Tanuma 
et al. 199£ ; Zimmer ? Lesch, &: Birk 1997J), and pro- 



impo'tant role as observed in solar flares by X-ray bata fc Yokoyama 1999 ). The magnetic reconnection 



satellite Yohkoh ([Tsuncta 1996; Shibata 1996|), and in 



substorms by satellite such as Geotail (e.g. Hoshinc 



occurs at numerous current sheets in Galactic disk 
and halo, which may be called "galactic flare" (|Stur-| 



et al. 199£). In the solar atmosphere, the magnetic rock &: Stern 198C|). Recently, Tanuma et al. (1999) 



reconnection heats the plasma from a temperature of 
several xlO 6 K to several xlO 7 K (sometimes sev- 
eral xlO 8 K), and accelerates it to 10 2 ~ 3 km s _1 (e.g. 
Shibata 199€ ; Yokoyama fc Shibata 1995J ) . If an ex- 
plosion such as a supernova compresses the current 
sheet, magnetic reconnection occurs and heats inter- 
stellar gas. 

Historically, the steady reconnection mechanisms 
were proposed on the basis of analytical studies. In 
the Sweet(1958)-Parker(1957) type reconnection, the 
diffusion region is so long as to occupy whole cur- 
rent system. Therefore, the stored magnetic energy 
is released by Ohmic heating, which can hardly be 
applied to solar flare phenomena, because the recon- 
nection rate of the Sweet-Parker model is too small 

— 1/2 

(~ R m ) in the solar corona, where R m is magnetic 
Reynolds number (~ 10 10-15 ). On the other hand, 
in the Petschek(1964) type reconnection, the diffusion 
region is localized near an X-point, and standing slow 
shocks occupy whole current systems. In the config- 
uration, the motor effect associated with slow shocks 
is much more dominant than Ohmic heating, which 
can hence be applicable to solar flare phenomena, be- 
cause the reconnection rate of the Petschek model is 
~ 0.1 — 0.01, independent of R m when R m is large. 
This is called "fast reconnection". Two models have 
been suggested to clarify how the fast reconnection 
can be realized. One is so-called "externally driven 
fast reconnection model" , which predicts that the fast 
reconnection should be controlled by boundary condi- 
tions (Petschek 1964; Sato &: Hayashi 1979J ; Priest & 
Forbe s 1986 ). The other is so-called "spontaneous 
fast reconnection model", which predicts that the 
fast reconnection mechanism spontaneously devel- 
ops from inside the system by self-consistent interac- 
tion between plasma microscopic processes (current- 
driven anomalous resistivities) and macroscopic re- 
connection flows (see also [Ugai 198*6]) . 

The purpose of this paper is to examine the mag- 
netic reconnection triggered by a supernova by per- 



performed numerical simulations and suggested that 
the magnetic reconnection can generate the Galactic- 
ridge X-ray emission (see section |4.2j ). In this paper, 
we solve the 2D resistive MHD equations numeri- 
cally to see how the magnetic reconnection occurs 
when a supernova triggers the magnetic reconnection 
in Galaxy. In this paper, we analize the detailed time 
evolution of the reconnecting current sheet, study de- 
pendence of the results on many parameters, and ap- 
ply the results to Galaxy. 

In the next section, we describe the method of the 
2D resistive MHD numerical simulations of the mag- 
netic reconnection triggered by a supernova as the 
simplest model. In section 3, we describe the results 
of the simulations, and we apply the results to the 
X-ray plasmas in Galaxy in section 4. In the last 
section, we summarize this paper. 



2. METHOD OF NUMERICAL SIMULATION 

2.1. The Situation of Problem 

The magnetic reconnection occurs when the inter- 
stellar magnetic field collides with another magnetic 
field which is not exactly parallel. The magnetic re- 
connection heats and accelerates interstellar gas by 
releasing magnetic energy. In the present paper, we 
treat the magnetic reconnection triggered by a super- 
nova (figure |l]) . We assume a simple initial condition 
(uniform temperature everywhere, and uniform mag- 
netic field outside the current-sheet) and large simu- 
lation region. We assume the magnetic fields are anti- 
parallel to each other, and a supernova occurs near 
the current sheet. We solved the two-dimensional 
(2D) resistive magnetohydrodynamic (MHD) equa- 
tions numerically to see how the magnetic reconnec- 
tion heats the interstellar gas. 



2.2. MHD Basic Equations 
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The 2D resistive MHD basic equations are written 
as follows: 



dp 

~dt 



+ V • (pv) = 0, 



(1) 



p— + p{v ■ V)v + Vp g = -J x B + pg, (2) 



dt 



dB 

~dt 
de 



V x (v x B) = -cV x (r]J), 



(3) 



0/ 



+ V- [(e+p g )w] = n\J\ 2 + v- Vp g (4) 



where /?, v, B, n, e, J, and g are mass density, veloc- 
ity, magnetic field, electric resistivity, internal energy, 
current density (= cV x B/4ir), and gravity (=0), re- 
spectively. We use the equation of state for the ideal 
gas, i.e., p g = (7 — l)e where 7 is the specific heat 
ratio (=5/3). 

2.3. Normalization 

We normalize length, velocity, and time. The units 
are 



H^WOElNPf 



-1/3 



pc, 



C s = (jR s T ) 1/2 ~ 20T 4 1/2 km s" 1 , 



(5) 
(6) 



-1/3 



T = g-3xl Si4 /3 T 4 - 1/2 (fj yr, (7) 

respectively, where E51, T4, P13, C s , and i? g are the 
initial explosion energy (-£7 ex ) in unit of 10 51 erg, tem- 
perature in unit of 10 4 K, gas pressure in units of 
10 -13 erg cm -3 , sound velocity, and gas constant, re- 
spectively. We assume that H is derived from 



E P 



Pc 



.(^Lrfi 



7-1 v 3 



H Z D SN 



(8) 



where D^n (= H) is a half-thickness of assumed ini- 
tial high-pressure region in a direction perpendicu- 
lar to x-y plane, and p cx = p Ta tioPg = 500p g = 
5 x Kr n pi3 (Model Al). 

The units of density, gas pressure, magnetic field 
strength, current density, and resistivity are po ~ 
1 x 1CT 25 g cm -3 , p C 2 ~ 4 x 10~ 13 erg cm" 3 , 
{poC 2 /^) 1 ' 2 ~ 3.2 pG, {pqC^/Stt^/H, and 
c 2 HC s /A-7r, respectively. In this paper, we solve non- 
dimensional basic equations, and describe the results 
of numerical simulations by non-dimensional values. 

2.4. Anomalous Resistivity Model 

We assume the anomalous resistivity model as fol- 
lows: 



'to 



% + a(vd/v c ~ 1) if v d > v. 



if v d < v c , . 

2 it „,. ^„, l y J 



([Ugai 1992| ; |Yokoyama fc Shibata 1995) ), where v d (= 
J/p), p, J, v c , and 770 are the normalized rela- 
tive ion-electron drift velocity, non-dimensional mass 
density, non-dimensional current density, threshold 
above which anomalous resistivity sets in, and "back- 
ground uniform resistivity". The parameters are 
rjQ = 0.015, a = 0.1 and v c = 100 in this paper for 
the typical model (Model Al). 

2.5. Initial Condition of Typical Model (Model Al) 

Figure [2] shows the initial condition of our nu- 
merical simulation. We assume the temperature is 
T = Tq = 1 (uniform) everywhere. The current-sheet 
thickness is 2/ mit = 2. We take Cartesian coordinate 
(x,y). The magnetic field, gas pressure, and density 
are 

B = B t&nh(y)x, (10) 

Pg = P0 + Pmag COsh~ 2 (y), (11) 

P = 1P%I T = Po + ilPm%/ T o) cosh~ 2 (y), 

respectively where Bq : 
0.6, po = 1, and p££* 
(= B 2 /8ir = po/P) outside the current sheet. The ra- 
tio of gas to magnetic pressure is j3 = 8ttpo/Bq = 0.2 
(\y\ > 1). Total pressure (P = p g + |S| 2 /87r) is 
uniform everywhere. The sound velocity is C s = 
(iPg/p) 1 = 1 (uniform) in the initial condition. The 
initial Alfven velocity is iff* = 5 /(4vrpo) 1/2 =* 2.45 

(lvl>i). 

We assume the symmetric boundaries for the top 
(y ~ 81.7) and bottom (y ~ —54.8) surfaces, and 
the periodic boundaries for the left (x ~ —120.9) and 
right (x ~ 120.9) surfaces. The simulation region size 
is L x ~ 241.8 times L y ~ 136.5. The number of grid 
points is (N x x Ay) = (465 x 602). The grid sizes are 
non-uniform (Ax > 0.20, Ay > 0.025). The number 
of grid points in the current sheet is 80. We put a 
high-pressure region instead of a point explosion (su- 
pernova) at (x,y)=(0, y ex )=(0, 7), whose radius and 
gas pressure are r$ = 1 and p cx = p ra ,tioPg = 500po> 
respectively. 

We use the 2-steps modified Lax-Wendroff method. 
We neglect gravitational force, effects of rotation such 
as shearing motion and Coriolis force, cosmic rays, 
viscosity, radiative cooling, and heat conduction in 
our simulation models. The cooling times due to heat 
conduction and radiation are 



12) 



8.68, x = (1,0), p = 1/7 = 

nj is initial magnetic pressure 



7"cond 
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X 



r rad 



v 100 pc 

nkT 
n 2 A(T) 

10 3 ( ' 



) 2 ( 



10 4 K' 



-5/2 



yr, 



(13) 



■)( 



ii 



10 4 K /v 0.1 cm" 3 ' 



A(10 4 K) 



n -l 



10 21 erg cm 3 s 1 



yr, 



(14) 



respectively, for cool, dense gas, where A is the 
cooling function ( Spitzer 1962| ), and kq is constant 



(= 10 6 erg s 1 cm x K i ). The radiative cooling 



time becomes r rac j ~ 10 6 (n/0.1 cm 3 ) 1 yr for the 
hot gas of T ~ 10 7 K, where A(10 7 K) is ~ 10~ 23 erg 
cm 3 s _1 . Furthermore, in the X-ray gas near Galac- 
tic plane (Galactic Ridge X-ray Emission=GRXE; 
T ~ 10 8 K and n ~ 3 x 10~ 3 cm -3 ; see Discus- 
sion), the radiative cooling time is ~ 10 10 ' 5 yr, which 
is much longer than the typical time scale (3 x 10 6 
yr; see next section) of the physical process exam- 
ined in this paper, so that the radiative cooling can 
be neglected. On the other hand, the conduction 
cooling time for the X-ray gas is ~ 30 yr [which is 
~ 10 5 (A cff /l kpc)(T/10 8 K)- 1 / 2 yr in fact because 
GRXE gas is collisionless] <3x 10 6 yr, so that the 
conduction cooling can not be neglected, where A c fr 
is the effective length of helical magnetic field. Nev- 
ertheless, we neglect the effect of heat conduction for 
simplicity in this paper, since the basic properties of 
magnetic reconnection such as reconnection rate nor 
energy-release rate are not much affected by the heat 
conduction ( Yokoyama fc Shibata 1997| ). We assume 
so long current sheet that the effect of evaporation 
is negligible, though the evaporation occurs if the in- 
terstellar cloud exists near the reconnection region 
and the reconnected field lines penetrate the cloud. 
In this paper, we assume the uniform strong field and 
long current sheet under pressure equilibrium for sim- 
plicity, instead of a small turbulence or small current 
sheets under non-equilibrium. 

3. RESULTS 

3.1. Typical Model (Model Al) 

Figure [| shows the time evolution of the two- 
dimensional distribution of the current density, den- 
sity, and gas pressure, with magnetic field lines, in an 
early phase of the interaction between the supernova 
and current sheet. A supernova occurs near the cur- 
rent sheet. A blast shock passes and perturbs the cur- 
rent sheet (figure ^) . Figures ||-^ show the time evolu- 
tion of spatial distributions of the temperature (figure 
H|), gas pressure (figure ||), and current density (fig- 
ure H), with velocity vectors and magnetic field lines. 



tion results. The current sheet evolves and magnetic 
reconnection occurs sequentially as the following 5 
phases, (i) Shock wave passage (t < 2; figure [7|a]) 
and the current sheet thinning by the tearing insta- 
bility in nonlinear phase (t ~ 100 — 200; figure [7][b] ) . 
(ii) Sweet-Parker type reconnection (t ~ 200 — 270; 
figure 0[c]). (iii) "Secondary tearing instability" in 
the Sweet-Parker current sheet (t ~ 250 — 270; figure 
0[d]). (iv) Generation of magnetic islands (t ~ 270; 
figure Qe]) and Petschek type (fast) reconnection 
with two standing slow shock regions (t ~ 270 — 400; 
figure 0[f]). (v) End of Petschek type reconnection 
(t~400). 

The current sheet is unstable to the tearing mode 
instability QFurth, Killeen, fc Rosenbluth 1963| ). The 
magnetic dissipation time, Alfven time, and time 
scale of the tearing instability in the initial state are 



init 
T dis 



init 2 



V 



T A 



init 



7? 0.015 

pnit 1 



67, 



r^j r^j 



2.5 



0.4, 



,_init /'^init^init N il/2 r o 



(15) 

(16) 

(17) 



respectively, where l imi (= 1) is the half-thickness 
of initial current sheet. The magnetic Reynolds 
(Lundquist) number is 



omit 
m,y 



.init /init 



H 1 



167. 



>/o 



(18) 



Phase I: Shock Wave Passage (t < 2) and 
Current-Sheet Thinning by Tearing Instabil- 
ity (t ~ 100 — 200) The supernova blast shock (a 
fast-mode magnetohydrodynamic shock) crosses the 
current sheet in an early phase (t < 2.5; figure ||). 
The magnetic reconnection does not occur immedi- 
ately after the shock wave crosses the current sheet, 
because the crossing time is too short to drive the 
magnetic reconnection. 

Figure || shows the time variation of half-thickness 
of the current sheet, defined by the minimum half- 
width of the half-maximum of gas pressure calculated 
parallel with y-axis in \x\ < 25. Figure ^ shows its 
schematic illustration. The magnetic field is dissi- 
pated (t ~ — 100), and the current-sheet thickness 
increases slowly. The current-sheet thinning occurs 
by the tearing instability (t ~ 100 — 200). Fig- 
ure [TTj shows the time variation of magnetic flux, 
defined by ip(t) = / * max^^o \(B x v y )[aA(x, 2)] + 
(B x v y )[at(x, — 3)]|dt. The process in this phase cor- 
responds to the nonlinear stage of the tearing-mode 
instability ([Furth. Killeen, fc Rosenbluth 1963] ; [Mag 



Figur |e~p] shows the schematic illustration of simula- ara fc Shibata 1999 ; Steinolfson and &: Hoven 1983 
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1984). The growth rate is consistent with the re- 
sults of Magara & Shibata (1999). The dashed line 
in figure [u] shows log^(t) ~ 0.00257*. It is 1/2 
times that of the linear phase [log ?/>(£) = log 10 e x 



0.43(l/r t init )(l/4vr)t ~ 0.0042/;; |Magara & Shibata 
1999[} It explains the result of our numerical sim- 
ulation well. 

The half-thickness of the current sheet decreases to 
I ~ 0.35 at t ~ 200 when the current-sheet thinning 
stops. The current-sheet length is comparable to the 
most unstable wavelength of the tearing instability 
occurring in the initial current sheet, i.e., 



A t , ~ 5.6iC" 1/4 r nit 



m,y 



(19) 



( |Magara k. Shibata 1999[ ). It is ~ 20 in this model. 
The current-sheet thickness in this phase is explained 
by 

k ~ ly)«m,t > 



(20) 



i.e., this current sheet corresponds to Sweet-Parker 
current sheet ( Magara &: Shibata 1999J ). Here, the 
magnetic Reynolds number is 



R 



m,t 



(At/2)< it 



m 



(21) 



The equations (|20|) and (|2~l|) give us l t ~ 0.23 and 
-Rm,t ~ 1467, respectively. These values explain our 
results in this phase well. In summary, we can say 
that the current-sheet thinning stops when the cur- 
rent sheet evolves to Sweet-Parker current sheet. The 
gas flows into the diffusion region at the velocity of 
V[ n ~ 0.01. The duration of this phase is At t ~ 100, 
which is approximately equal to the time for ini- 
tial current sheet to become the Sweet-Parker sheet, 
l init /v in ~ 1/0.01 ~ 100 ~ 20r t init . 

Phase II: Sweet-Parker Type Reconnection 



(t ~ pOO- 270) 0weet(195S)-Parker(1957) type re- 



conn eclion slaiLs aL L ><> 200 in an asyninieliic silua- 
tion. Alfven velocity is higher in the region {y > 0) 
because of the low density gas generated by the su- 
pernova. The gas flows along current sheet, and is 
accelerated to v ~ 2.5, which is equal to the Alfven 
velocity near (x, y) ~ (0, —2). Figure 11 shows (a) the 
time variation of velocity of the gas flowing into the 
reconnection region. The inflow velocity is defined by 
v m = max| a .| <10 [v y [at(x, -3)] - v y [at(x, 1)]] /2. Fig- 
ure [y](b) shows Alfven velocities attained at (x$,2) 
(the solid line) and at (xo,— 2) (the dashed line) 
where xo is the point where maximum inflow veloc- 
ity attained. Figure |i~l|(c) shows the reconnection 
rate defined by v- m /vA,s, where va,s is the smaller 
one between VA[at(a^,l)] and vp\at(x$, — 3)]. The 



current-sheet length increases to A ~ 80 (figure ||), 
because the gas in the initial current sheet flows out 
along the current sheet. The inflow velocity becomes 
v- m ~ 0.025. The ratio of inflow to outflow velocity 
(~ 0.025/2.5 ~ 0.01) is explained by 






Asp/2 



R 



-1/2 
m,SP 



(22) 



where i? m ,SP is the local magnetic Reynolds number 
defined by 

(Xsp/2)v a 



R 



m,SP 



'/0 



(23) 



Asp and lsp are the length and half-thickness of the 
Sweet-Parker current sheet ( Parker 1957 ), respec- 
tively, at the time when the secondary tearing insta- 
bility starts . These equations give us v\ n /v A ~ 0.01 
with -R mi sp ~ 10 4 , Asp ~ 80, and /sp ~ 0.3, respec- 
tively, in this model. The condition for the setup of 
the secondary tearing instability is determined by the 
aspect ratio of Sweet-Parker current sheet, so that the 
secondary tearing instability occurs when aspect ra- 
tio [l ength/thicknes s= (A/2)//] becomes larger than 
100 ( |Biskamp 1993 ). The equation ( f22|) explains the 
result of our simulation in this phase {t ~ 220 — 250) 
well. During the Sweet-Parker type reconnection, 
the current-sheet thickness increases in proportion to 
A 1 ' 2 , which is derived from the equation (|22|). 

At t ~ 200 — 220, the inflow velocity (vm) increases 
to ~ 0.15 (see the arrow in figure |ll| ), which is much 

— 1/2 

larger than Sweet-Parker value ~ R S p v A ~ 0.025. 
This transient inflow velocity is determined by f; n ~ 
2/2io w out/At ~ 0.15, which is derived from the mass 
conservation equation A t Wi n ~ 2/210^0111) assuming the 
density is constant and w out (~ 0.6) is the velocity 
of the outflowing gas along the current sheet at this 
time (t ~ 210), where Z210 ~ 2.5 is the current-sheet 
half-thickness in x ~ 20 (at the edge of Sweet-Parker 
curre nt sheet) at t ~ 210 (Cf. |Fu fc Lee 1986| ; |Shibata 
T997| ). 



Figure ^ shows the time variation of various ener- 
gies in the reconnection region (\x\ < 25 and \y\ < 3). 
The magnetic energy (the solid line), thermal en- 
ergy (the dotted line) , and kinetic energy (the dashed 
line) are defined by E m3ig = J_ 25 f_ 3 (.B 2 /87r)dxdy, 



Eth = J-25/-3bg/(7 ~ l)]dxdy, and E kin = 
J_25 J-3(/ 9t ' 2 /2)dxdy, respectively. They are cor- 
rected by the energy flux at the boundaries (|x| = 25 
and \y\ = 3). Figure [l3| shows the schematic illus- 
tration of the time variation of the magnetic energy- 
release rate. The magnetic energy is converted to the 
thermal energy slowly by the tearing instability and 
Sweet-Parker type reconnection. At t ~ 200 — 220, 
the magnetic energy-release rate increases (see the ar- 
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rows in figures |l2] and 13) by the nonlinear growth of 
tearing instability. 

Phase III: Secondary Tearing Instability (t ~ 
250 — 270) The Alfven time and time scale of the 
tearing instability of the Sweet-Parker current sheet 
become, 



7"A,SP 



Tt,SP 



^-0.12, 

%^A 



0.85, 



(24) 
(25) 



respectively. The tearing instability (we call "sec- 
ondary tearing instability" in this paper) occurs in 
Sweet-Parker current sheet at t ~ 250. The current- 
sheet thickness decreases again. Some magnetic is- 
lands are generated with intervals of ~ 5 — 10 in x, 
which are determined by the wavelength of secondary 
tearing instability, 



A 2 -4.9</ 4 y / S p, 



(26) 



where -R m ,y is the local Reynolds number, defined by 



R 



vaL 



SP 



m,y 



m 



(27) 



and ua is the local Alfven velocity. The equation (26) 
gives us A2 ~ 4 with R m ,y ~ 50, which explains the re- 
sults of our simulation in this phase well. Some small 
islands collide with one another, and are coalesced to 
become big islands. The secondary tearing instability 
goes on for At2t ~ 20. It is equal to 2.5r t sp ~ 20. 



while it is ~ 0.02 during Sweet-Parker type reconnec- 
tion. In this simulation, the energy-release rate dur- 
ing Petschek type reconnection is only a few times 
higher than that during Sweet-Parker type reconnec- 
tion. The latter, however, is much small in the ac- 
tual Galaxy because it is determined by the magnetic 

Reynolds number as e = v- in /v & = R m . The gas in 
the diffusion region are accelerated along the current 
sheet by magnetic tension force to v ~ 2.5 — 3, which 
is comparable to the Alfven velocity in the region 
where the maximum velocity of reconnection jet is 
attained. 

The magnetic energy is converted to the ther- 
mal energy slowly during the tearing instability and 
Sweet-Parker type reconnection, and quickly during 
Petschek type reconnection, particularly while mag- 
netic islands are moving along the current sheet from 
the reconnection region (figure [IJ). The energy con- 
version is suppressed when the magnetic islands are 
staying in the diffusion region (t ~ 310 and 360). 
Total released magnetic energy is ~ 5500 in the re- 
connection region (|x[ < 25 and \y\ < 3; figure 12), 
which is much higher than the initial thermal en- 
ergy of the supernova (E e ^ ~ 1440). The maximum 
energy release rate and maximum heating rate are 
-dE mag /dt ~ dE th /dt ~ 70 - 80 at t ~ 340. They 
are determined by the Poynting flux entering the re- 
connection region, 



dE ma g d 



(29) 



Phase IV: Petschek Type Reconnection (t ~ 
270 — 400) Figure |lj shows the time variations of 
(a) maximum drift velocity (vd = J/p) in the cur- 
rent sheet (|xj < 5 and —2.6 < z < 1.0), (b) density 
where maximum v^ is attained, (c) current density 
(J) where maximum v^ is attained, and (d) maxi- 
mum magnetic reconnection rate defined by r/J. The 
density decreases from 0.5 to 0.1 at the center of the 
current sheet (figure |lj) because the gas flows away 
along the current sheet from the dissipation region. 
The maximum drift velocity exceeds the threshold 
(v c = 100) in the reconnection region immediately af- 
ter the big magnetic island moves away. The anoma- 
lous resistivity sets in at t ~ 270. The half-thickness 
of the current sheet decreases to I ~ 0.15 (figure ||). 
The diffusion region is much localized near the X- 
point. Petschek (1964) type reconnection begins with 
various violent phenomena. The gas flows into recon- 
nection region at v m ~ 0.1 — 0.2 (figure |ll|). The 
magnetic reconnection rate (figure O) increases to 



— ~ 0.05 - 0.10, 

VA 



(28) 



where S is the reconnection region size. The above 
equation gives us ~ 80, which explains the results of 
the simulation (The reconnection region size is de- 
rived as S ~ 37 if we use the values of B ~ 8 and 
v m ~ O.2.). 

The current density increases in the current sheet 
(at t = 210 and 250; figure |6|) and behind the moving 
magnetic islands (at t = 275 and 300; figure ^). The 
slow shocks are formed behind the moving islands. 
Two standing slow shocks are generated in the cur- 
rent sheet (at t = 325 and 350; figure ||). Figure 15 
shows the profile of the various physical quantities in 
x = — 12 at t = 325. The gas pressure, temperature, 
density, and velocity are high in the current sheet. 
The magnetic field is weak in the current sheet. The 
current density is large at the slow shock regions so 
that the profile of current density has double peaks. 
It is noted that this Petschek-type reconnection is 
very time-dependent; Some magnetic islands are re- 
peatedly generated and ejected in the current sheet, 
and magnetic energy-release rate is controlled by the 
magnetic islands. 
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Phase V: End of Fast (Petschek Type) Re- 
connection (t ~ 400) The current-sheet length is 
limited by the high gas pressure near the right and left 
symmetric boundaries (x = —120.8 and x = 120.8). 
Figure [l6| shows the distribution of gas pressure in 
whole simulation region. The current sheet, whose 
half-thickness was I = 1 at t = 0, expands near 
the left and right boundaries. The magnetic loops are 
accumulated at the left and right boundaries so that 
current-sheet thickness increases from ~ 2 at t = 160 
(figure [L6|[a]) to ~ 20 at t = 300 (figure Mb]), and 
eventually to ~ 50 at t = 450 (figure |l§c]). The 
gas pressure increases in the reconnected loop. The 
current-sheet length decreases to A ~ 50 at t ~ 450 
(figure Il6|[c]). The ejection from the diffusion region 
is suppressed by the high gas pressure in the loop near 
the boundaries, so that the global field configuration 
tends to be the potential field produced by the local- 
ized currents near the boundaries. On the other hand, 
the magnetic tension force suppresses the inflow to 
reconnection region because the potential magnetic 
field outside the current sheet is bent. Consequently, 
the magnetic reconnection rate decreases suddenly at 
t ~ 380. Duration of the Petschek type reconnection 
is At-p ~ 170. 

3.2. Parameter Survey 

3.2.1. Dependence of Results on Magnetic Field 
Strength ((3) 

We examine the dependence of results on the mag- 
netic field strength [B Q = (8irp /(3) 1/2 ], i.e., (3 (Model 
A2 and Bl-5). The time variations of magnetic 
energy-release rates are shown in figure |l7| by the 
dotted line (Model Bl; j3 = 0.1), solid line (Model 
A2; 13 = 0.2), and dashed line (Model B2; (3 = 0.3). 
The magnetic energy-release rate decreases with (3 as 
~ 250 (p = 0.1), ~ 70 {(3 = 0.2), and ~ 20 {(3 = 0.3), 
respectively. 

Figure [l8| (a) shows the /3-dependence of maximum 
magnetic energy release rates during Petschek type 
reconnection (the solid line) and during Sweet-Parker 
type reconnection (the dotted-line) in the reconnec- 
tion region (|x| < 25 and \y\ < 3). The energy-release 
rate decreases with (3. The rate during Petschek type 
reconnection is determined by the Poynting flux en- 
tering into reconnection region, 



"-C'mag I r,D Q a— 3/2 

— tHMp ~ 2— Sv m OC (3 *', 

dt 47T 



(30) 



where S is the reconnection region size, v- m = cva = 
eB/(4:TTp) 1 ' 2 is the inflow velocity to the reconnection 
region, and e is the reconnection rate. This theoreti- 
cal /3-dependence (30) is also shown in the figure |l~8|(a) 
by the dashed line [~ 2( J B 2 /4vr)Se< it ~ 6.8/T 3 / 2 



where S is derived to be ~ 26 if we use the results of 
the typical model]. It explains the numerical results 
well. On the other hand, the magnetic energy-release 
rate during Sweet-Parker reconnection is given by 



dE, 



mag 



dt 



|SP ~ 2^-Sv m oc /T 5 / 4 , (31) 

47T 



where S is the reconnection region size, v- in = 

— 1/2 1/2 

R m gpt>A oc v A is the inflow velocity to the recon- 
nection region. This theoretical /^-dependence (J3l|) is 
also shown in the figure |l8|(a) by the dashed-dotted 

line [~ 2(B$/4Tr)SR^ 1 Hvf it ~ 2.5/T 5 / 4 where S is 
derived to be ~ 67 if we use the results of the typical 
model]. It explains the numerical results well. 

Figure |l^(b) presents the times when maximum 
magnetic energy release rates are attained during 
Petschek (the solid line) and Sweet-Parker (the 
dashed line) type reconnections. The time when 
Sweet-Parker type reconnection starts is determined 
by the time scale of dissipation in an early phase 
(~ Ai t ) plus the time scale of tearing instability 
(~ At t ). It is almost equal to ~ 2At t - The time scale 
of tearing instability (r t au A oc /J 1 ' 4 ) is presented 
in figure |l8|(b) by the dashed-dotted line (~ 230/3 1 / 4 ). 
The time when Petschek type reconnection starts is 
determined by the time when Sweet-Parker type re- 
connection starts (~ 2Ai t ) plus the time for the 
Sweet-Parker current sheet to get long until it be- 
comes unstable to the secondary tearing instability, 
A/2 > 100/sp ~ 30 ( |Biskamp 19931 ), which is equal to 
~ Asp/va- The Alfven time (ta oc v^ oc /3 1 ' 2 ) is pre- 
sented in figure |l8|(b) by the dashed line (~ 580/3 1 ' 2 ). 

Figures |i~8|(c) and (d) show the maximum temper- 
ature of heated gas and maximum velocity of the re- 
connection jet (\y\ < 2). The maximum temperature 
and the maximum velocity decreases with (3. The in- 
terstellar gas is heated, by passing the reconnection 
region, to 



T, 



max - I 1 + S I \ J ■'■ m 
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(32) 



which is derived from the equation n out kT max ~ 
ninkTi n + B 2 /8ir. Here, we assume (ni n /n out )Ti n is 
independent of (3. The /3-dependence is plotted in 
figure |l8|(c) by the dashed line [~ 7(1 + 1/(3)]. In 



the simulation, a part of hot gas of the supernova 
flows into the reconnection region, and re-heated by 
the reconnection. The gas is accelerated to Alfven 
velocity, 



B 



"jet ~ VA 



(4tt^) 1 /2 



oc/3" 



-1/2 



(33) 



by the magnetic tension force along the current sheet, 
which is shown in the figure |i~8|(d) by the dashed line 
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(~ 1.5/3 1 ' 2 which is 1.5 times higher than u™ lt be- 
cause va increases with time). 

3.2.2. Dependence of Results on Distance Between 
Supernova and Current Sheet 

We examine the dependence of results on the dis- 
tance (y ex ) between the supernova and current sheet 
(Models A2 and Cl-3). Figures |l|(a) and ||(b) il- 
lustrate the maximum magnetic energy release rate 
and the time when it is attained, respectively. The 
magnetic energy-release rate decreases with y cx . The 
time scale increases with y cx . It is because the initial 
perturbation decreases with y ex , the time when mag- 
netic reconnection starts increases and more magnetic 
field dissipates before the reconnection starts. The 
maximum temperature (figure |l9|[c]) decreases with 
y ex . The maximum velocity of reconnection jet (fig- 
ure 19 [d]) does not depend on the distance because 



it is determined by the Alfven velocity. 

3.2.3. Dependence of Results on Initial Energy of 

Explosion 

We put the point explosion with the gas pressure 
of E ex oc p ex = p ra tioPo a s the initial condition. We 
examine the dependence of results on p ra tio oc E ex 
(Models A2 and Dl-5; figure |2Q) , The normaliza- 
tions given by the equations (|5])-(|7j) change if p ra tio of 
the equation (0) changes. We must rewrite the units 



of length and time as H oc p rat ( and r oc p vat { , re- 



1/3 



-1/3 



spectively. Figure |2C](a) presents the maximum mag- 
netic energy release rate. It does not depend on the 
energy of initial explosion. Figure 2C(b) shows the 
time when maximum magnetic energy release rate is 
attained. It decreases with p ra tio because the initial 
perturbation decreases. Figure |20|(c) shows the max- 
imum temperature of gas heated by the reconnection. 
The temperature increases with p ra tio because the hot 
gas of the initial supernova is re-heated by the mag- 
netic reconnection. Figure |20|(d) shows the velocity of 
reconnection jet. It does not depend on p ra tio because 
it is determined by the Alfven velocity. 

We examine the model in which two explosions 
appear symmetrically in both sides [at (0, 7) and 
(0, -7)] of the current sheet (p ratio = 500; Model D5). 
The results are plotted in figure [l8| by the cross ( x ) . 
The time scale is much shorter because the initial 
current sheet is compressed from both sides. Other 
results are similar to those of typical model. 

3.2.4. Dependence of Results on Resistivity Model 

We examine the dependence of the results on the 
background uniform resistivity (tjq) and the resistiv- 
ity models (either anomalous or uniform for whole 
evolutions). Figure [2l] shows the time variations of 



the half-thickness of the current sheet, for Model A2 
(t/o = 0.015; thick solid line), Model Fl (r/ = 0.03; 
thin dotted line), Model F2 (rj = 0.0075; thin 
dashed-dotted line), Model F3 (r/ = 0.00375; thin 
dashed-dotted-dotted line), Model F4 (r? = 0.001; 
thin dashed line), Model F5 (rj = 0.0001; thin solid 
line), and Model F6 (77 = 0; thick dashed line). In the 
current sheet of Model A2 (770 = 0.015), the tearing 
instability {t ~ 60 — 150) is followed by Sweet-Parker 
type reconnection {t ~ 150 — 200), secondary tear- 
ing instability (t ~ 200 — 240), and Petschek type 
reconnection (t ~ 240). In zero-resistivity model 
(Model F6), the current sheet does not dissipate at 
first. It, however, starts to become thin at t ~ 220, 
and eventually leads to Petschek type reconnection 
by numerical resistivity (resistivity-like effect by nu- 
merical noise). The result of Model F5 (770 = 0.0001) 
is similar to that of Model F6 (77 = 0), not that of 
F4 (r/o = 0.001). The numerical resistivity in these 
simulations is estimated as r? num = 0.001 — 0.0001. 
In Model Fl (77 = 0.03), the current sheet dissipates 
(t < 100), and the magnetic reconnection does not 
occur (t < 800). The maximum temperature of the 
gas heated by the reconnection and the velocity of 
the reconnection jet do not depend on tjq. 

The result of uniform resistivity model (Model E) 
is also shown in FigurdTq by triangle (A). Petschek 
type reconnection does not occur in the uniform re- 
sistivity model. The magnetic energy-release rate 
and the time scale are approximately equal to the re- 
sults of anomalous resistivity model during the Sweet- 
Parker type reconnection. 

3.2.5. Dependence of Results of Simulation Region 
Size and Grid Size 

We examine the dependence of results on the simu- 
lation region size (L y ) in y axis (Models A2 and Gl- 
2). The magnetic energy release rate increases with 
the simulation region size because the stored initial 
total magnetic energy increases. The other results do 
not depend on the simulation region size. 

We examine the effects of grid size (Ax) (Models 
Al-2). Until the secondary tearing instability starts, 
the time evolution of current sheet depends on the 
parameters such as (3 and r\ well so that the time 
scales of phenomena are predictable. At this stage 
(the onset of the secondary tearing instability) , how- 
ever, it becomes difficult to predict where small is- 
lands appear and when and how they move and col- 
lide with one another, since the current-sheet thick- 
ness decreases to a few x Ay after the secondary tear- 
ing instability starts. At this stage, the current sheet 
is too thin to be free from numerical resistivity in our 
problem. The maximum energy release rate, maxi- 
mum temperature of the heated gas and velocity of 
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the reconnection jet do not depend on Ax, although 
the time when maximum magnetic energy release rate 
is attained depends on Ax (figure 22). 



4. DISCUSSION 

Odstrcil & Karlicky (1997) performed an interest- 
ing 2D numerical simulation of the magnetic recon- 
nection triggered by a point explosion near the cur- 
rent sheet, with parameters different form ours, in 
order to construct a model of the solar flare trig- 



gered] by a disLanL flaie. — They assumed iesislivil i y 
as i] oc | J — J c |, where J c is Hie critical currenL den- 
sity above which the anomalous resistivity sets in. 
Zimmer et al. (1997) tried to explain that the gas 
is heated in Galactic halo by magnetic reconnection 
triggered by the collision of the high-velocity clouds 
(HVC; M ~ 1O 7 - 8 M ; V ~ 100 km a -1 ; [Blitz 1999| ), 
by performing the 2D MHD numerical simulation. 
Birk et al. (1998) performed 2D numerical simulation 
of the magnetic reconnection including the ionization 
and recombination. Above three groups followed the 
current-sheet thinning and the generation of the mag- 
netic island. They, however, did not reveal the time 
evolution of the secondary tearing instability, because 
they did not use enough small grid size. In this pa- 



per, ive use 80 grid puinls, su LhaL we levealed huw 
the 



secondary tearing inslabilily leads Lu Pelsehek 
type magnetic reconnection. 

We revealed how the fast reconnection starts at the 
current sheet via the secondary tearing instability by 
assuming the anomalous resistivity model (equation 
||) in this paper. Main result of this paper is re- 
vealing how the current sheet leads to fast reconnec- 
tion via the secondary tearing instability. This result 
does not depend on the form of resistivity before the 
anomalous resistivity sets in ( |Yokoyama &: Shibata| 
1995 ; Ugai 1992 ). The simulation results we apply to 
the actual Galaxy, such as temperature of heated gas, 
do not depend on the resistivity model ( |Yokoyama &| 
Shibata 1995] ; |Ugai 1992Q . 

We assume that the initial current-sheet has en- 
hanced gas pressure [p g = Po + (po//?)cosh~ (y); 
B x = Z?otanh(y); B y = 0] in the paper. The tear- 
ing instability grows similarly in both this model 
(enhanced gas pressure model) and force-free model 
bg = Po; B x = B cos0.5iT(y + 1); B z = 
J Bosin0.5vr(y + l); B y = 0; |Magara fc Shibata I999J . 
The basic physics of magnetic reconnection do not 
change, though the nonlinear phase starts earlier in 



the enhanced gas pressure models ( Magara &: Shi 



bata 1999| ) . The physical processes of our results are 
hardly influenced by the boundary conditions because 
the simulation region is large enough. We do not as- 
sume the initially enhanced resistivity in the current 
sheet to see how the spontaneous reconnection oc- 



curs after the secondary tearing instability. We use 
so small grid size that we can suppress the numerical 
resistivity in our numerical simulations. 

The results do not depend on p ra tio very much. 
The magnetic reconnection occurs with any p ra tio- 
It means that the magnetic reconnection occurs not 
only by a supernova, but also by the small per- 
turbation. The magnetic reconnection can be trig- 
gered by various mechanisms such as a supernova 
(or a point explosion), superbubbles (sum of su- 
pernovae; Tomisaka 1998; [rcnorio- Tangle fc Bodci~ 
hcimer 1988J ), stellar winds, collision of interstellar 



clouds, shock waves, Parker(1966) instability (undu 
lar mode of magnetic buoyancy instability), magne- 
torotational instability ( palbus &: Hawley 199 1| ), col- 
lision of galaxies, cosmic rays etc. in numerous cur- 
rent sheets. 

The random field stores the magnetic energy at ~ 
several 10 pc ( Rand fc Kulkarni 1989J ). Total mag- 
netic energy is £ mag ~ [{B)' 2 ohs /8Tr] V G ~ 10 544 erg, 
where (-B) bs ~ 3 pG is assumed as the mean ob- 
served field strength, and Vq ~ 2 x 10 2 kpc 3 is the 
volume of Galaxy. The thermal energies of GRXE 
(Galactic Ridge X-ray Emission) and coronal halo 
are -Egrxe r -' 10 55 erg (e.g. [Kaneda et al. 1997|) 



and E halo ~ 10 572 erg ([Habe fc Ikcuchi 198q ; |Li& 



Ikeuchi 1992; Ikeuchi 1988| ), respectively. Total rota- 
tional energy of Galaxy (-E^ot ~ 10 58 - 9 erg) and total 
kinetic energy of the interstellar gas (-Eism ~ 10 58 ' 2 
erg) are higher than the thermal and magnetic ener- 
gies ( [Sturrock Sz Stern 1980| ). The rotational energy 
can be converted to the thermal energy by Galactic 
dynamo and magnetic reconnection in ~ 10 8 yr (e.g. 
iParker 197l] ; fTanuma et al. 1999|) . 



4.1. Implication for Observations 

In our numerical simulations, the magnetic 
Reynolds number is much smaller than actual 
value in Galaxy. The Alfven velocity, resistiv- 
ity, and magnetic Reynolds number are va ~ 
6 x 10 6 (B/W /iG)(n/0.1 cm" 3 )" 1 ^ cm s -i j v „ 
10 7 (T/10 4 K)- 3 / 2 cm 2 s -1 , and R m ~ Xv A /ri ~2x 
10 20 (5/10 /xG)(n/0.1 cm- 3 )- 1 /2(A/l00 pc)(T/10 4 K) 3 / 2 , 
if we assume the length of current sheet which is get- 
ting thin by the tearing instability, where the local 
magnetic field strength, number density, and temper- 
ature are A ~ 100 pc, B ~ 10 pG, n ~ 0.1 cm -3 , and 
T ~ 10 4 K, respectively. Sweet-Parker current-sheet 
thickness is 



lsp ~ Ait! m 1/2 
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The time scale of Sweet-Parker reconnection is 
A A 

T"SP 



^in V A R m 1/2 

~3x io i6 (^^) 3 / 2 (— — r 1/2 

V 100pc ; V 10//G ; 

(7Tr^-^)- 3/4 (T7^7) 3/4 yr (35) 



"0.1 cm 



TO 4 K' 



(under classical Spitzer resistivity), so that the Sweet- 
Parker type reconnection can not play an important 
role as energy-release mechanism in actual interstel- 
lar medium. 

On the other hand, the time scales of the secondary 
tearing instability and current-sheet thinning are 



r 2 , tear ~(— ) i/2 ~ 4) 3/2 (-x 2 , 

Vva A v A 
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respectively. The secondary tearing instability oc- 
curs faster than the current-sheet thinning if T2 : tear < 
Tthiru i-e., the current-sheet thickness decreases to 

' < ^2,tear 

= XR- 1/3 
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In this situation, the aspect ratio becomes 
A/(2Ztear) ~ 4 x 10 6 . Hence, the secondary tearing 
instability starts before Sweet-Parker current sheet 
develops completely. The time scale and wavelength 
of the secondary tearing instability are calculated to 
be 



T 2 ,tear ~ 2 X 10 6 ( 
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respectively. 

When a plasmoid moves along the current sheet, 
the inflow velocity increases behind it so that the 
current-sheet width decreases to 
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(41) 



by assuming A 2 ,tear'Vm = k,teaxVA, where 77 in is the in- 
flow velocity toward the dissipation region. In this 
transient large inflow stage, the anomalous resistiv- 
ity can set in, which eventually leads to Petschek 
type reconnection, since the condition of anoma- 
lous resistivity is that the current-sheet thickness 
must be smaller than the ion Lamor radius [~ 
10 7 (T/10 4 K) 1 / 2 (S/10 fiG)- 1 cm; figure ||. The 
current-sheet thickness, otherwise, decreases to the 
ion Lamor radius through third, forth, and fifth tear- 
ing instabilities (figure 23). The third tearing in- 



stability occurs when the current-sheet thickness be- 
comes ~ 1 x 10 12 cm (wavelength is ~ 7 x 10 15 cm; the 
time scale of the third tearing instability is ~ 10 yr) 
in the current sheet generated by the secondary tear- 
ing instability. The forth, fifth, sixth, seventh, and 
eighth tearing instabilities start when the current- 
sheet thickness becomes ~ 4 x 10 11 cm, ~2x 10 10 
cm, ~ 1 x 10 9 cm, ~ 7 x 10 7 cm, and ~4x 10 6 cm, 
respectively. Total released magnetic energy is 



32 D 

\3 . m507 D \2, 



A 



£~(f)A 3 ~10 50 (-f-) 2 (— — ) 3 erg, (42) 
8-7T 10 /jG 100 pc 

if we assume the volume of space filled with the mag- 



netic 



assume tne volume ol space niled witn tne mag- 
field is A 3 . Petschek reconnection continues for 

A 



rp 



ev A 



A 



which is equal to one-third of the time of Galactic 
rotation. The magnetic energy-release rate is 

^^l 3 4.5 ( ^ r l ( ^^_ )2 

r P v 0.1 7 v 100pc' 

10 fiG 0.1 cm d 
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4.2. Generation of X-ray Gas in Galaxy 

Figure |24| shows the 2D distribution of the thermal 
radiation [F oc n 2 A(T)], for the typical model (Model 
Al), where A(T) is the cooling function (~ T a ). The 
above a is 0.5, —0.6, 2.9, and 1.5, respectively when 
the unit of temperature is To ~ 10 7 K (figure |j[a]), 
~ 10 5 ~ 6 K (figure |§b]), 10 4 K (figure ||[c] ) , and 10 3 
K (figure ^[d]). The X-ray thermal emission shown 
in figures p4|(a) and 24(b) may be observed by Newton 
and Chandra in the galaxies. 

The bright thermal X-ray emission are observed 
along the Galactic plane (Lx ~ 10 38 erg s -1 ; Galac 



tic Ridge X-ray Emission=GRXE; e.g. Kaneda et 
al. 19971; |Koyama et al. 1986a|; |Makishima 1994|; [Ya 



masaki et al. 1997; Yamauchi et al. 1996; Valinia & 



Marshall 1998). If the magnetic reconnection oc- 
curs in the magnetic field of B ~ 30 /iG, which 
is generated locally by Galactic dynamo, the hot 
(T ~ 10 8 K) component (n ~ 3 x 10 -3 cm~ 3 ; 
|Kaneda et al. 1997 ) of GRXE can be generated 
( [Makishima 1994 |Tanuma et al. 1999|) . After then, 
the hot gas is confined by the locally strong mag- 
netic field for r cond ~ 10 5 (A cff /l kpc)(T/10 8 K)" 1 / 2 
yr, where A e ff is the effective length of helical 
magnetic field. Note that the radiative cooling 
time is T rad ~ 10 10 - 5 (T/10 8 K)(n/0.003 cm" 3 )- 1 
[A(10 8 K)/10~ 23 erg cm 3 s -1 ] -1 yr, which is too 
much longer than r cont j so that the radiation is neg- 
ligible, where A(T) is the cooling function. The 
heating time by Petschek type reconnection is rp ~ 
5 x 10 5 (A/100 pc)(£/30 fiG)' 1 ^/ 0.003 cm" 3 ) 1 / 2 yr, 
so that the heating is balanced with the conduction 
cooling, where A is the thickness of the reconnecting 
magnetic field. 

The hot gas are also observed in Galactic halo 
flPietz et al. 199q ), Galactic center (e.g. Koyamaj 
|1989| ), and clusters of galaxies (e.g. |Sarazin 1986| ). 
The magnetic reconnection can generate the hot 
gas ( Zimmer, Lesch, & Birk 1997; Birk, Lcsch, & 



Ncukirch 1998|; |Kerp, Lesch, fe Mack 1994| ; |Kerp et 



irp 

m 



al. 1996; Soker 1997|) and confine it in the magnetic 



loop for a long time. 

The reconnection jet is accelerated in opposite di- 
rections by the magnetic tension force, which is like 
bipolar-jet. Inside the helical tube and current sheet, 
the temperature and gas pressure increases by the 
self-pinch of helical magnetic tube ot the collision be- 
tween the reconnection jet and interstellar cloud so 
that the star formation may be triggered. 

5. SUMMARY 

We examined the magnetic reconnection triggered 
by a supernova (or a point explosion) in interstel- 
lar medium, by performing two-dimensional resis- 
tive magnetohydrodynamic (MHD) numerical simu- 



lations with high resolution and large simulation re- 
gion size. We found that the magnetic reconnection 
starts long after a supernova shock (fast-mode MHD 
shock) passes a current sheet. The long current sheet 
evolves as follows: (i) The supernova shock perturbs 
the current sheet. The magnetic reconnection does 
not occur immediately after the passage of the su- 
pernova shock across the current sheet, because the 
magnetic field lines do not have enough time to recon- 
nect. The tearing-mode instability is excited by the 
supernova shock, and the current sheet becomes thin 
in its nonlinear stage, (ii) The current-sheet thin- 
ning is saturated when the current-sheet thickness 
becomes comparable to that of Sweet-Parker current 
sheet. After that, Sweet-Parker type reconnection 
starts, (hi) The current-sheet length increases during 
the Sweet-Parker type reconnection, because the gas 
in the Sweet-Parker current sheet flows out along the 
current sheet. The secondary tearing mode instabil- 
ity occurs in the Sweet-Parker current sheet, because 
the ratio of length to thickness of the current sheet 
exceeds ~ 100. (iv) Further current-sheet thinning 
occurs after the islands moved away. The anoma- 
lous resistivity sets in, and Petschek type reconnec- 
tion starts. The magnetic energy is released while 
the magnetic islands are moving during the Petschek 
type reconnection. The phenomena depend only on 
/? (i.e. B), not on the distance between the supernova 
and current sheet, nor the energy of initial explosion. 
In actual Galaxy, the magnetic reconnection can 
occur when the non-parallel magnetic fields (i.e. cur- 
rent sheet) exist and a point explosion such as su- 
pernova, superbubble, collision of the clouds or stel- 
lar wind will occur near the current sheet. The ini- 
tial current sheet will be made by the turbulence, 
Parker(1966) instability ( [Tanunia et al. 1998]), mag - 
netorotational instability QBalbus fc Hawley 1991 ), 
collision of galaxies etc. The magnetic islands (or he- 
lical tubes) confine the hot gas, and flow fast along 
the current sheet. The interstellar gas flows also in- 
side the magnetic helical tube in three-dimensional 
(3D) view. The gas pressure is very high inside the 
magnetic islands so that the gas flows along the tube 
actually. In the magnetic tube and current sheet, 
the temperature and gas pressure are very high so 
that the star forming region may appear. We will ap- 
ply the results to the X-rays from the Galactic ridge, 
Galactic halo, Galactic center and clusters of galaxies. 
We will examine the reconnection in 3D simulations 
in future works. 



The authors thank K. Makishima in University of 
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Fig. 1. — The situation of our problem. The interstellar space is filled with the magnetic fields. The magnetic fields are not 
exactly parallel to one another, and have numerous current sheets at all scales. We examine how the magnetic reconnection occurs 
when a supernova perturbs the current sheet. We calculate the interaction between the supernova and current sheet by assuming 
simulation region shown by black square. 

Fig. 2. — Initial condition of the typical model (Model Al) of our numerical simulations. We assume the uniform temperature 
and uniform total pressure (P g + B 2 /8iv) everywhere. The current-sheet thickness is 2P nlt H = 2H where H is the unif of 
length. The units of velocity and time are sound velocity (Cs) and H/Cs, respectively. We solve the non-dimensional equations 
numerically in this paper. The magnetic field is assumed as B x — Bo tanh(y//f) and B y = 0. The plasma (3 is 8tvpo/Bq = 0.2 
everywhere outside the current sheet. The distance between the supernova and current sheet is y cyi H — 7H. Minimum grid size 
is (Ax, Ay) > (0.20#, 0.0255). The number of grid points is (465 x 602). The simulation region is -120. 9H < x < 120.9// and 
— 54.8// < y < 81.7H. We assume no gravity (g — 0). 



Fig. 3. — Time evolution of 2D distribution of the current density (J), density (p), and gas pressure (p g ), with the magnetic 
field lines, in an early phase for the typical model (Model Al). The half-thickness of initial current sheet is I — 1 (\y\ < 1). 
The initial density, gas pressure, and magnetic field outside the current sheet are no — 1, po — 0.6, and Bq ~ 8.7, respectively. 
The plasma j3 is 8-rrpo/B 2 — 0.2 (\y\ > 1). The Alfven velocity is i? /(47rp)" 1/2 ~ 2.5. We put a supernova at (0, y ex ) = (0,7). 
The supernova expands in the magnetic field. These figures show the expanding supernova shock well. The slow shock propagates 
along the magnetic field. The shock front is in a high-density region, while the center of supernova is in a low-density region. 
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Fig. 4. Time evolution of 2D distribution of the temperature (T), magnetic field lines, and velocity vectors, in a main phase 
of the magnetic reconnection triggered by a supernova, for the same model shown in figure ti (Model Al). A supernova appears 
at (x,y) = (0, 7), and disturbs the current sheet. The tearing instability occurs, and the current sheet gets thin (t ~ 100 — 200). 
Sweet-Parker type reconnection occurs at t ~ 200 — 250. Petschek type reconnection starts at t ~ 270. The magnetic reconnection 
heats the gas, and accelerates it by the magnetic tension force. The hot gas of supernova flows into the diffusion region. The 
magnetic islands are generated in the current sheet sequentially. The heated gas is confined in the magnetic island. 



Fig. 5. — Time evolution of 2D distribution of the gas pressure (p g ), magnetic field lines, and velocity vectors, in a main phase 
of the magnetic reconnection triggered by a supernova, for the same model shown in figure [J (Model Al). The initial gas pressure 
is p g = po + (po//3) cosh - (y) where po = 0.6 and (3 = 0.2. The gas pressure is high in the current sheet, especially inside magnetic 
islands (i ~ 270, 300). 



Fig. 6. — Time evolution of 2D distribution of the current density (J), magnetic field lines, and velocity vectors, in a main phase 
of the magnetic reconnection triggered by a supernova, for the same model shown in figure (Model Al). The initial magnetic 
field is B x = Bo tanh(y) and B y — where Bo = 8.58. The plasma /3 is po/(Bg/87r) = 0.2. The initial current-sheet thickness is 
2l lnlt = 2. The initial current density is ~ Bo/2 ~ 4.5 in the current sheet. The current density increases in the diffusion region. 
Petschek type reconnection generates two standing slow shocks in the current sheet (see also figure FL3) . The slow shocks appear 
also behind the moving magnetic islands (see the figures at t = 270 and 300). 



Fig. 7. — The schematic illustration of the results of our simulation (Model Al). A supernova disturbs the current sheet (figure 
a). The tearing instability occurs long after the passage of fast shock wave (figure b). The current-sheet thickness decreases. 
Sweet-Parker type reconnection occurs (figure c). The secondary tearing instability occurs in the Sweet-Parker current sheet 
(figure d). A magnetic island is generated (figure e). Petschek type reconnection starts after the islands moves away (figure f). It 
heats the interstellar gas, and accelerates the heated gas by the magnetic tension force to the Alfven velocity. Many islands are 
generated sequentially, and confine the heated gas. 



Fig. 8. — Time variation of the half-thickness of the current sheet (I), which is defined by the minimum half-width of half- 
maximum of the gas pressure calculated parallel with y-axis in \x\ < 25, for the same model shown in figure |2| (Model Al). The 
supernova shock wave perturbs the current sheet (t < 2). The magnetic field dissipates slowly (t < 100). The tearing instability 
occurs (t ~ 100 — 200). The current-sheet thinning occurs during the tearing instability. Sweet-Parker type reconnection occurs 
(t ~ 200 — 250). The secondary tearing instability occurs in the Sweet-Parker current sheet (t ~ 250 — 270). The current-sheet 
thickness decreases again by the secondary tearing instability. Petschek type reconnection occurs (t ~ 270 — 400). The current 
sheet gets much thin during the Petschek type reconnection. At t ~ 390 — 400, the current sheet is very thin because the magnetic 
island is moving. 
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Fig. 9. — The schematic illustration of the time variation of the current-sheet thickness. At first, the magnetic field dissipates 
slowly. Second, the current-sheet thinning occurs by the tearing instability in a nonlinear phase. Sweet-Parker type reconnection 
occurs in the thin current sheet. Next, the secondary tearing instability occurs in the Sweet-Parker current sheet. The current- 
sheet thickness decreases more. Petschek type reconnection starts. The solid line moves to upper side for large- /3, small- j/ ox , and 
large-p ox models, while it moves to lower side for small- /3, large- y ex , and small-p ex models. 



Fig. 10. — Time variation of the magnetic flux, for the same model shown in figure H (Model Al). The magnetic flux is defined 
by ^{t) = Jo max H<io \(B x v y )[at(x, 2)] + (B x v y )[at(x, — 3)]|dt. The magnetic flux increases during the tearing instability. The 
magnet ic flux grows by log ib (t) ~ 0.00257£, which is shown by the dashed line. The rate is 1/2 times the growth rate of the linear 



phase (Magara fc Shibata 1999) 



Fig. 11. — Time variation of the inflow velocity toward the diffusion region, for the same model shown in figure 3 (Model 
Al). Figure (a) shows the inflow velocity, defined by Vin = maxj^io [uy[at(a;, —3)] — v y [at(x, 2)]] /2. Figure (b) shows the Alfven 
velocities at (xo, 2) (dotted line) and (:rn, —3) (solid line), where xo is the point where maximum inflow velocity is attained. Figure 
(c) shows the reconnection rate defined by e = «i n /«A,s, where «a, s is smaller one between i>A[at(:ro, —2)] and VA[at(xo, 2)]. The 
inflow velocity and reconnection rate are Vi n ~ 0.1 — 0.2 and e ~ 0.05 — 0.1 during Petschek type reconnection, while they are 
~ 0.05 and ~ 0.02 during Sweet-Parker type reconnection. At the beginning of Sweet-Parker type reconnection (t ~ 200 — 220), 
they increases to ~ 0.16 and ~ 0.6, respectively, because of the nonlinear growth of the tearing instability (see the arrow). 



Fig. 12. — Time variation of the various energies, showing how the energy conversion occurs by the magnetic reconnection in 
a central region (|a;| < 25 and |y| < 3) for the same model shown in figure (Model Al). The magnetic energy is converted to 
the thermal energy slowly by Sweet-Parker type reconnection and secondary tearing instability (t ~ 200 — 270), and quickly by 
Petschek magnetic reconnection (t ~ 270 — 360). The magnetic islands suppress the magnetic reconnection when they are staying 
at the reconnection region (t ~ 310 and 360). The magnetic energy-release rate increases to ~ 40 at the beginning of Sweet-Parker 
reconnection because the inflow velocity increases (see the arrow). Maximum energy release rate is — d_B mag /di ~ 70 — 80 during 
Petschek type reconnection (t ~ 340). The non-dimensional thermal energy of the initial supernova is E ex ~ 1800 in the numerical 
simulation. Main thermal energy is supplyed from the magnetic energy, not from the initial supernova. 



Fig. 13. — The schematic illustration of the time variation of the magnetic energy-release rate in the reconnection region. 
The magnetic dissipation occurs slowly at first. The tearing instability starts much after the passage of supernova shock. Sweet- 
Parker type reconnection occurs after current-sheet thinning by the nonlinear phase of the tearing instability. The magnetic 
energy-release rate increases once at t ~ 200 — 220 (see the arrow) because the inflow velocity increases. After then, the magnetic 
energy is released slowly by Sweet-Parker reconnection. The secondary tearing instability occurs in Sweet-Parker current sheet. 
The magnetic energy-release is suppressed by the small islands generated by the secondary tearing instability. Petschek type 
reconnection starts at last. The magnetic energy is converted to the thermal energy violently during Petschek type reconnection. 
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Fig. 14. — Time variations of (a) the maximum drift velocity (vd — J/p), (b) density p where the maximum ^d is attained, (c) 
current density (J) where the maximum v c \ is attained, and (d) maximum reconnection rate denned by rjj (—2.6 < y < 1, \x\ < 10), 
for the same model shown in figure ti (Model Al). The current density increases to J ~ B/lsp ~ 20 after Sweet-Parker reconnection 
starts. The density in the current sheet decreases during Sweet-Parker type reconnection and secondary tearing instability 
{t ~ 200 — 270). The drift velocity exceeds the critical value (v c = 100) so that the anomalous resistivity sets in (t ~ 270). 
Petschek type reconnection starts so that the reconnection rate increases. 



Fig. 15. — Profile of the variables in x = —12 at t ~ 342, for the same model shown in figure H (Model Al). p g , T, p, and 
|w y | are higher in the current sheet than those outside the sheet. The magnetic field strength is lower in the current sheet. The 
profile of the current density has two peaks (y ~ —2.3 and 0.7), corresponding to the slow shocks. The density ratio across the 
shock front is X = niln\ ~ 0.3, where subscripts 1 and 2 mean that the values are attained at y ~ —2.1 and —0.3, respectively. 
The as pressure ratio is p g 2/p g i ~ 0.1, while 0.47 is required by Rankine-Hugoniot relation. 



Fig. 16. — 2D distribution of the gas pressure at t — 160 (a), 300 (b), and 450 (c), for the same model shown in figure H (Model 
Al). Figures (a), (b), and (c) show the distributions during the tearing instability, a little after Petschek type reconnection starts, 
and after Petschek type reconnection stops. The gas is ejected from the diffusion region toward the left and right boundaries. The 
gas pressure increases in the current sheet near the boundaries. It extends the current sheet near the boundaries. Furthermore, 
the magnetic tension suppresses the magnetic reconnection because the magnetic field is bent outside the current sheet. 



Fig. 17. — Time variation of the magnetic energy-release rate (\x\ < 25, \y\ < 3|), for Model Bl (/3 = 0.1; dotted line), Model 
A2 (/3 = 0.2; solid line), and Model B2 (/3 = 0.3; dashed line). The magnetic energy-release rate decreases with /3. The maximum 
magnetic energy release rates are — d_E mag /di ~ 250, 70, and 20, respectively. 



Fig. 18. — Dependence on /3 of (a) the maximum energy release rate (d_B mag /di) in the reconnection region (\x\ < 25, \y\ < 3), 
(b) the time when maximum d_E mag /dt is attained, (c) maximum temperature of heated gas, and (d) maximum velocity of the 
reconnection jet (Model A2, Bl-5, and E). (a) The magnetic energy-release rate decreases with j3. The solid line and dotted 
line show the maximum magnetic energy release rates by Petschek and Sweet-parker type reconnections, respectively. The rates 
are determined by Poynting flux entering into the reconnection region. The theoretical /3-dependence is shown by the dashed 
line [2(_Bo/47r)5'euX llt °c f3~ 3 where we assume e ~ 0.1 and S ~ 26] for Petschek type reconnection, and the dashed-dotted line 
[2(Bo/4i")Sii~ gp^A 1 ' of fl~ 5 where we assume i? m ,sp ~ 10 4 and S ~ 67] for Sweet-Parker type reconnection. The results of 
the uniform-resistivity model (Model E) are also shown by the triangle (A), (b) The time scale increases with /3. The magnetic 
phenomena occur in the Alfven time (ta oc /3 1 , dashed line) or the time scale of tearing instability (ta oc /3 1 , dashed-dotted 
line), (c) The maximum temperature decreases with j3. It is determined by T max oc (1 + 1//3). (d) The maximum velocity decreases 
with (3. The gas is accelerated to Alfven velocity va oc t^ 1 oc f3~ 1 ^ 2 . 
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Fig. 19. — Dependence of the results on the distance (j/ex) between the supernova and current sheet (Models A2 and Cl-3; 
see the caption of figure Ko). Figure (a) shows the maximum magnetic energy release rate in the reconnection region (\x\ < 25, 
\y\ < 3|). It decreases with y cx . Figure (b) shows the time when maximum magnetic energy release rate is attained. It increases 
with y ox because the initial perturbation decreases. Figures (c) and (d) show the maximum temperature of the heated gas and 
maximum velocity of the reconnection jet. They do not depend on y cx . 



Fig. 20. — Dependence of the results on the energy of the initial explosion (E ex oc p ex ~ Pratio) (Models A2 and Dl-4; see the 
caption of Figure [IS). Figure (a) shows the maximum magnetic energy release rate in the reconnection region (jx| < 25, \y\ < 3). 
It does not depend on p pr atio. Figure (b) shows the time when maximum magnetic energy release rate is attained. It does not 
depend on p ra tio. Figure (c) shows the maximum temperature of the gas heated by the magnetic reconnection. It increases with 
Pratio because the hot gas flows into the reconnection region. Figure (d) shows the velocity of the reconnection jet. It does not 
depend on p pra tio because the velocity of the reconnection jet is determined by Alfven velocity. The results of two-supernovae[at 
(0, 7) and (0, — 7)]-model (Model D5) are indicated by x. The magnetic reconnection starts early in the model because the current 
sheet is compressed from both sides. 



Fig. 21. — Time variations of the half-thickness of the current sheet (/), defined by the minimum half-width of half-maximum 
of the gas pressure in \x\ < 25, for Model A2 (770 = 0.015; thick solid line), Model Fl (770 = 0.03; thin dotted line), Model F2 
(770 = 0.0075; thin dashed-dotted line), Model F3 (rjo = 0.00375; thin dashed and double-dotted line), Model F4 (r/o = 0.001; thin 
dashed line), Model F5 (r? = 0.0001; thin solid line), and Model F6 (r? = 0; thick dashed line). In Model A2 (770 = 0.015), the 
tearing instability (t ~ 60 — 150) is followed by Sweet-Parker type reconnection (t ~ 150 — 200), secondary tearing instability 
(t ~ 200 — 240), and Petschek type reconnection (t > 240) (figure H). In zero-resistivity model (F6), the current sheet does not 
dissipates. It starts to become thin at t ~ 220, and Petschek type reconnection occurs by only numerical resistivity (resistivity-like 
effect by the numerical noise). In Model Fl (rjo = 0.03), the current sheet dissipates early (t < 100). 



Fig. 22. — Time variations of the magnetic energy-release rate in the reconnection region (\x\ < 25, \y\ < 3), for Model Al 
(Aa; = 0.20; solid line) and Model A2 (Aa; = 0.25; dotted line). The magnetic energy-release rate does not depend on the grid 
size because the "background uniform resistivity" (770) is enough large to suppress the numerical resistivity. The time scale of 
phenomena changes because the grid number of the initial supernova changes, and because the perturbation on the current sheet 
changes. 



Fig. 23. — The schematic illustration of the time variation of the current-sheet thickness in actual interstellar medium. Here we 
assume B ~ 10 /iG, n ~ 0.1 cm" 3 , T ~ W 4 K, and Spitzer conductivity. The current-sheet thickness decreases to the ion Lamor 
radius (~ 10 7 cm) by the nonlinear phase of secondary tearing instability or by successive tearing instabilities. When it decreases 
to the radius, Petschek type reconnection can occur in the interstellar medium. 



Fig. 24. — 2D distribution of the thermal radiation (F oc n 2 T a ) for the same model shown in figure H (Model Al). a is 0.5, 
—0.6, 2.9, or 1.5, respectively, when the unit of temperature is assumed to be To ~ 10 7 K (figure [a]), ~ 10 5 " 6 K (figure [b]), 10 4 
K (figure [c]), or 10 3 K (figure [d]), The high-temperature gas in the island and current sheet is bright. 
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Table 1 

Parameters. We studied the effect of the grid size (Models Al-2), plasma /3 (Models Bl-5), distance y ox between the current sheet and 
supernova (Models Cl-3), energy of initial supernova (Models Dl-5), resistivity model (Models E and Fl-4), and simulation region size 

(Models Gl-2). Typical model is Model Al. 



Model 



Wcx 



Pratii 



'/o ■ 



SN' 



Ax% 



Ay h 



Lx 1 



N x 



Ny l 



Al 


0.20 


7.0 


500 
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a the ratio of the gas to magnetic pressure outside the current sheet 
b the distance between the supernova and current sheet 
c the initial gas pressure of the supernova 
d "background uniform resistivity" 

e the parameter of anomalous resistivity; 0.1 for the anomalous resistivity, and 0.0 for the uniform 

resistivity 

' number of the initial supcrnovae 
s minimum grid size in :r-axis 
h minimum grid size in y-axis 
'simulation region size in x-axis 
J simulation region size in j/-axis 
k number of grid points in ir-axis 
'number of grid points in j/-axis 
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